Changing shapes in the nanoworld 
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What are the mechanisms leading to the shape relaxation of three dimensional crystallites ? Kinetic 
Monte Carlo simulations of fee clusters show that the usual theories of equilibration, via atomic 
surface diffusion driven by curvature, are verified only at high temperatures. Below the roughening 
temperature, the relaxation is much slower, kinetics being governed by the nucleation of a critical 
germ on a facet. We show that the energy barrier for this step linearly increases with the size of the 
crystallite, leading to an exponential dependence of the relaxation time. 



Imagine a world where marbles would fuse upon con- 
tact, just as two water droplets usually do to minimize 
their surface energy and reach their equilibrium configu- 
ration. This peculiar behavior is thought to be usual in 
the nanoworld, because experiments suggest that objects 
in the. nanometer range can change shape in reasonable 
timesErcl, even if they are solid. This fact is crucial for the 
production and control of nanostructures, since these are 
generally obtained in out-of-equilibrium conditions, and 
are therefore metastable. The rapid shape relaxation of 
these solid particles (which lack the collective atomic dif- 
fusion mechanisms found in the liquid state) is dominated 
by surface diffusion of their surface atoms. Atomic diffu- 
sion is random (brownian) but nevertheless generates a 
global mass transfer from the high curvature regions (of 
higher chemical potential, roughly because atoms have 
less neighbors there) to the low curvature regions. For 
very small objects, surface diffusion is very efficient as a 
mass transfer mechanism. However, it has to be pointed 
out that this whole picture assumes that a continuous 
description of the objects is valid, thus allowing the defi- 
nition of a curvature-dependent local chemical potential. 
This is correct as long as the particle is large enough (to 
allow defining a chemical potential) and has a disordered 
or rough surface (in order that the chemical potential 
is diffcrentiable) as for a liquid droplet. .Within these 
assumptions, Herring, Nichols and MullinsEl have shown 
that this mass transfer mechanism leads to an equilib- 
rium time t eq which increases as the fourth power of the 
object linear size, which helps understanding why macro- 
scopic objects are kinetically allowed to violate thermo- 
dynamics with total impunity. It is important to note 
that this fourth power law is extensively used to predict 
equilibrium times fits, to extract diffusion constants from 
equilibration ratesou. In this paper, we address the im- 
portant question of what happens when the temperature 
is below the roughening transitionB, so that the particle 
is faceted, a common situation for nanostructures. There 
is no general agreement on a continous approach in this 
regimeQ, and therefore we use Monte Carlo simulations 
to show that the equilibration time of a faceted parti- 



cle does not simply vary as predicted above. Indeed, 
we show that, although surface diffusion remains the im- 
portant matter transport channel, the equilibration time 
is limited by another physical mechanism - facet nucle- 
ation. This different physical mechanism leads to a very 
different size dependence of t eq , limiting the validity of 
the continuous law and explaining why nanobjects might 
spend much more time in metastable equilibrium than 
expected. 



Physical model 



We use standard Kinetic Monte Carlo simulations!] to 
study the equilibration of unsupported 3D crystallites 
having a perfect fee crystalline structure. Since we are 
only interested in finding generic laws for the size de- 
pendence of t eq (which should not depend on the details 
of atom-atom interaction), we have chosen a very simple 
energy landscape for atomic motiora. We assume that 
the potential energy E p of an atom is proportional to its 
number i of neighbors, and that the kinetic barrier E act 
for diffusion is also proportional to the number of initial 
neighbors, before the jump, regardless of the final number 
of neighbors, after the jumptl : E act — —E p — i*E where 
E sets the energy scale (E — 0.1 eV throughout the. pa- 
per). Comparing with recent ab-initio calculationscj for 
the Al(lll) surface, we note that our one-barrier assump- 
tion does give the good order of magnitude of the rela- 
tive jump frequencies for the different hopping processes 
of interest here. We. also exclude any explicit "Ehrlich- 
Schwoebel" barriertil for atoms hopping around corners. 
Therefore, the probability pi per unit time that an atom 
with i neighbors moves is pi = vq exp[— i * E/ki,T], where 
iso = 10 13 s^ 1 is the Debye frequency. Thus, our sim- 
ple kinetic model includes only one parameter, the ratio 
E/ksT where fee is the Boltzmann constant and T the 
absolute temperature. We have changed the tempera- 
ture from 300K to 800K, and the crystallite number of 
atoms from 700 up to 13000. The initial configuration 
of the clusters is elongated (same initial aspect ratiola) , 
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and we stop the relaxation when the crystallites are close 
to equilibrium, with an aspect ratio of 1.2. 



Simulation results 



Fig. [I] shows a log-log plot of the relaxation time as a 
function of the number of atoms in the crystallite. The 
continuous law predicts a slope of 4/3, which agrees with 
our simulations only for the highest temperature, 800AT. 
As the temperature decreases, the slope continuously in- 
creases, reaching much higher values than this 4/3. This 
strong deviation from the continuous law suggests that 
an altogether physical mechanism limits the mass trans- 
fer at low temperature. The constant increase of the 
exponent is a clue that an exponential dependence of t eq 
on size and temperature might be present. 

To investigate the different behaviors at high and low 
temperatures, we examine (Fig. ^|) the different mor- 
phologies of the clusters at T = 700A and T = 300A'. 
At high temperatures, many kinks and steps can be seen, 
indicating that the continuous approximation for the cur- 
vature might be valid. On the contrary, at low tempera- 
tures the crystallite is fully faceted, with angular points 
and edges, making it difficult to define a chemical poten- 
tial properly. Moreover, the presence of facets makes it 
impossible to transfer atoms from the cluster tips to its 
central region by simple atomic diffusion : atoms reach- 
ing the facets do not find a trapping site there and even- 
tually get back to kinks or steps in the tip regions. This 
means that the crystallites can be trapped for long times 
in these faceted configurations at low temperatures, as 
can be clearly seen on Fig. || : the cluster progressive ap- 
proach of equilibrium is continuous at high temperature 
and more steplike at low temperature, indicating that 
at high temperature atoms continuously attach to the 
central regions, whereas some more discontinuous mech- 
anism operates at low temperatures. A careful exami- 
nation of low temperature relaxation pictures suggested 
that the transition from one step (i.e. a fully faceted, rel- 
atively stable configuration) to the following (lower) step 
demands the nucleation of a germ on a large facet. This 
germ grows and eventually forms a new atomic layer, thus 
bringing the crystallite closer to equilibrium (see Fig. |^). 
Therefore, we expect that the limiting process for matter 
transfer at low temperatures is the formation of a critical 
nucleus, as suggested from classical nucleation theorytj. 

This physical picture can be.-aartly quantified using 
the kinetic theory of nucleationlii The time needed to 
build the critical nucleus is given by : 

AG* 

tnucl OC exp(— — ) (1) 

k B T 

where AG* is the free energy barrier the system has to 
cross. We stress that here the situation is more complex 
than in the gas-liquid transition, where the atoms form- 
ing the incipient liquid critical cluster all come from the 
gas phase, which has a fixed chemical potential. Here, the 



tip atoms come from different environments with differ- 
ent energies (see the different colors in Fig. |^), making it 
difficult to calculate an average chemical potential. The 
key point is however to examine whether AG* depends 
on the particle size, thus creating an exponential contri- 
bution to the size dependence of t eq . .— . 

We use a classical umbrella sampling techniquai-3 to 
compute the crystallite free energy as a function of the 
number of atoms in the nucleating germ. The umbrella 
technique consists in adding a bias potential to the hamil- 
tonian of the system to force it to stay in a configuration 
of interest, even if it is unprobable, as is the case here for 
nucleation of the germ- Fig. || shows that AG* increases 
for larger crystallites^, which implies (Eq. |l|) that the 
nucleation time (and therefore t eq ) depends exponentially 
on the size of the cluster (provided, of course, that the 
AG* increase is not logarithmic, see below). 

What are the microscopic mechanisms leading to this 
AG* increase with crystallite size? The free energy of a 
nucleating germ is given bylij : 

AG = 2j line ^q - qA/j. (2) 

where q is the number of atoms in the germ, 7;i ne the 
line tension of the germ, and Afi the chemical potential 
difference for an atom going from the tip to the facet. 
We fit the curves of Fig. || by Eq. || which gives ju n e 
and Afi. We find ju ne — 0.129 ± 0.013eV/atom, inde- 
pendent of the crystallite size. This value is, as expected, 
close to the binding energy. To understand the size de- 
pendence of AG* (which comes from Afi), one can, as 
a first approximation, treat the tips in a continuous way 
: assimilating them to half an ellipsoid, we can estimate 
the tip curvature. This gives a rough measure of the kink 
and step density on the tips, and therefore of the den- 
sity of more or less mobile atoms, which can contribute 
to mass transfer. With this approximation, and taking 
arbitrarily the atom chemical potential to be zero on the 
facet, we get for the chemical potential difference for an 
atom going from the tip to the facet : 

Afi = ^surface^ (3) 

where ^surface is the average surface tension on the tip 
and k its curvature. Finally, we obtain the free energy 
barrier for nucleation : 

2 2 i 

Afi 'Jsurface ^ 

Figure || shows that Eqs. || and |] are in good agree- 
ment with our simulations and give coherent values for 
Isurfacei close to the binding energy E — O.leV (Eq. 
I gives 0.179 ± 0.008 eV atom" 2 and Eq. | leads to 
0.175 ±0.035 eV atom' 2 ). 



Conclusion, Perspectives 



The physical picture of nanocrystallite equilibration is 
the following : above the roughening temperature, the 
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continuous approach works well and leads to the classic 
fourth power law, the mass transfer being via atomic dif- 
fusion from kinks or steps from the high curvature regions 
to the existing kinks or steps of the low curvature region 
(which act as traps). Below this temperature however, 
large facets do appear in the low curvature regions and 
no kinks or steps are available, preventing the diffusing 
atoms from sticking there. Therefore, the route to equi- 
librium has to involve nucleation of new atomic planes, 
which is much more difficult and needs more time, lead- 
ing to an exponential increase of t nuc (and therefore t eqi 
which is directly related) as a function of the cluster size. 
We actually see no reason why this picture would not ap- 
ply to much larger particles, up to the micrometer range. 
If defect-free particles of this size could be produced, be- 
low the roughening temperature they should show per- 
fectly flat facets which would demand the nucleation of 
a germ for effective mass transfer, thus generating an 
exponential-type size dependence, and probably prevent- 
ing any experimental observation of the equilibrationlij! 
Many open questions remain : the temperature depen- 
dence of t eq has to be understood, a more quantitative 
theory for AG* has to be worked out, and simulations 
in other geometries (including also the substrate) would 
also be of interest. 

We acknowledge useful discussions with H. Larralde, 
J.L. Barrat, L. Bocquet, J.J. Metois and J. Krug. 



atoms with more than 6 neighbors become very rare (there 
are always some particles with 6 neighbors which move 
much faster). 
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FIG. 1. Log- log dependence of the relaxation time as a 
function of the size of the crystallites for different tempera- 
tures. The slope of each linear fit is indicated. 

FIG. 2. Morphologies of crystallites of 1728 atoms at two 
different temperatures : (a) partially rough at 700 K and (b) 
fully faceted at 300 K. The colour of each atom depends on 
its number of neighbors. 

FIG. 3. Evolution of the total energy of crystallites as a 
function of the time logarithm for a 1728-atom cluster at two 
temperatures. The arrows in the low temperature curve in- 
dicate the transitions from one faceted configuration to the 
next. The total energy is defined as the number of atomic 
bonds times the bond energy (E=0.1 eV). At the end of each 
curve, the crystallite has almost reached its equilibrium shape. 
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FIG. 4. Cluster free energy during the formation of a nucle- 
ation germ on a facet as a function of the number q of atoms 
in the germ. The curves have been obtained at 400 K, for sev- 
eral cluster sizes which have approximately the same shape, 
close to equilibrium (their aspect ratio is indicated in the fig- 
ure). Clearly, the free energy barrier for the nucleation of a 
critical germ becomes larger as the crystallite size increases. 
Each solid curve is fitted by Eq. ^, allowing to obtain ju n e 
and A/i. 

FIG. 5. Dependence of the nucleation barrier (given by the 
maximum of the curves in Fig. ^|) on the tip curvature k 
(calculated from the crystallite shape). In the inset, we show 
the dependence of A^i (deduced from the fits in Fig. |^) on 
the curvature k. We fit these two curves by Eqs. ^| and |^ 
respectively, obtaining jsurface (see text). 



4 



-1 



D" 





& -3 
o 



-5 




• 800 K (1 .3) 
■ 700 K (1 .5) 

♦ 600 K (2.0) 
a 500 K (2.7) 

* 450 K (4.2) 
▼ 400 K (5.4) 



-7 



2.5 



3.0 



3.5 4.0 
log(N) 



4.5 



5.0 



